samp_n <- 100

rand_samp_data <- read.csv(file.path(dir_M, "git-annex/impact_acceleration/stressors/sst/int/rand_sample.csv"))



for(i in 1:samp_n){ # i = 2
 
plot_sst <- ggplot(dplyr::filter(rand_samp_data, sample_id==i), aes(y = sst_anom, x = year)) +
    geom_point(size = 2) +
   geom_line() + 
   stat_smooth(method=lm, se=FALSE, color="red")
  
plot(plot_sst)  
}